function y0=Hermite_interp(X,~,~,x0)
%X，Y是已知插值点向量序列
%DY是插值点处的导数值
%x0是插值点横坐标
%y0是待求的分段三次Hemie插值多项式在x0处的值% N表示向量长度
N=length(X);
for i=1:N
if x0>=X(i) &&x0<=X(i+1)
k=i;break;
end
end
al=x0-X(k+1);
a2=x0-X(k);
a3=X(k)-X(k+1);
y0=(al/a3)^2(1-2*a2/a3)*Y(k)+(-a2/a3)~2*(1+2*al/a3)*...
)+(al/a3)~2a2*DY(k)+(-a2/a3)~2*al*DY(k+1);